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Abstract 



This paper is an introductory guide to many-particle hydrodynamic interactions. Basic concepts 
of the fluid mechanics are assumed to be known. Experience in the Stokes equations is useful but not 
necessary. The study is estimated to fit five sessions about three hours each. 



Auguste Rodin: "Nothing else that I have done satis- 
fies me as much, because nothing else cost me so much 
effort..." Rj 
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1 Introduction 



Augustc Rodin: "...human thought is limited by compari- 
son with what nature transmits directly to us and imposes 
on us. All that is nesessary is to follow the model; char- 
acter results from its unity. . . "ry 



The text consists of three parts: a brief formulation of the problem, work sheets for own studies 
aiming to develop the basic concepts and a short concluding overview to indicate how these concepts are 
useful in construction of the theory modelling quantitatively many-particle systems [Feldcrhof 1976a]- 
| Cichocki 1995| . 



The structure of the curriculum emerged from the method of effective learning by inquiry [ McDcrmott ] , 
redesigned and extended according to our previous experience, to specific needs of the subject and to feed- 
back received from participants of a course guided by a preliminary version of the worksheets presented 
here. [ NfcDcrmott | , a non-standard textbook for teachers, called "a set of laboratory-based modules" by 
the authors, provides a method to develop habbits of effective learning, based on active inquiry, applica- 
tion of scientific reasoning and cooperative work in small groups. Connection with reality is essential in 
this learning pattern, since it provides the motivation leading to a personal engagement (students start 
from their own obsevations) and it serves as the natural objective evaluation of own understanding (stu- 
dents make experiments, which verify predictions of physical models which they have just constructed). 
A discrepency between own predictions and reality becomes a driving force to learn. An extension of this 
approach to study a theoretical cirriculum has been neither straightforward nor simplistic. Actually, it 
required a general analysis of creative learning principles. 

The course has been guided by a continuous struggle to take care of integrity of the learning process, 
namely to make its mission and its vision clear, its specific goals apparent, its structure simple and 
evolving to fit students' needs and capacities. Therefore we have started with a very specific formulation 
of the goal. We have tried not only to identify and to keep the right sequence of steps building on 
each other, but also to make this sequence apparent in advance. We have shown how to make various 
side connections, giving a chance to see possible generalizations and applications, to establish a relation 
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with own experience, and to recognize an own direction of further studies. The nesessary attitude of 
the instructor reaching out for integrity was to be first of all a student, challenged to develop a new 
deeper insight into learning/ teaching techniques as well as into physical and mathematical aspects of the 
hydrodynamic interactions, to make the own learning evident, and to allow other learners for influencing 
what and how they learn. This curriculum would be never developed without its practical application. 

These principles have led to the structure of the learning process presented here. The problem and the 
goal were formulated specifically in the announcement about the course (see Appendix) sent to scientists 
and students working on problems related to hydrodynamic interactions. This information was important 
to decide for participation in a non-standard activity. Originally the curriculum had been designed to 
be based only on active group work rather than passive listening to lectures. However, the participating 
scientists demanded to be conscious not only of a direction and goals, but also of a perspective of their 
studies, important in choosing what to investigate further. Therefore the structure has been modified. 
Work sheets (Sec. §-fh^) served as a guide in own studies carried out in small groups of 2-4 people during 



4 sessions about 2.5 hours each. An overview concluding lecture (Sec. 5.1-5.3) was added at the end 



of the course as a closure and as an application of the participants' own inquires. A similar pattern of 



education had been earlier developed and tested in [ Glazek,Maslowski, Wigckowski . 



Part I 

Formulation of the problem 

How to determine the behavior of N spheres in low Reynolds number incompressible fluid flow (N up to 
several hundred)? 

We will concentrate on the following 'friction problem': 
If translational and rotational velocities of the spheres are given, as well as an ambient fluid flow in which 
they have been immersed, then what are the forces and torques they exert on the fluid? 

This approach can be afterwards adjusted to solve also the twin 'mobility problem': 
If an ambient fluid flow and external fo rces and torque s acting on the spheres are given, then what are their 
translational and rotational velocities? Fclderhof 1988| 



Our goal is to inquire the basic structure and tools of the technique developed in | Felderhof 1976a 
Cichocki 1995| . 



Part II 

Developing basic concepts: work sheets 



Auguste Rodin: "I forced myself to express in each 
swelling of the torso or of the limbs the efflorescence of a 
muscle or of a bone which lay deep beneth the skin. And 
so the truth of my figures, instead of being merely super- 
ficial, seems to blossom from within to outside, like life 
itself. "Qj 



Following McDcrmott |, in this part we used different type styles to distinguish between a text guiding 
independent work (written like this sentence), general informations (slanted) and additional remarkes 
(small letters). 
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2 Principles of work 



The idea is to make the whole problem a subject of your own active inquiry, carried out and discussed in 
small groups of 2-4 persons, on the basis of work sheets written specially for you. 

Questions raised during our sessions will help to identify separate steps to be made, building subse- 
quently on each other. Each step consists of the problems (formulated as a separate subsection of the 
work sheets) to be solved by you. You may find it useful to keep a written record of your work. 

The end of each subsection is a point to conclude - first to share your reasoning with each other, and 
next to discuss your results with me, giving me your comments and questions. To keep track of time we 
will indicate in a 'calendar of progress' when your group has finished each subsection. 
If you have a problem blocking your progress, and any of your group cannot solve it, please ask me for 
help. 

Our interests and background vary. Therefore it is reasonable to divide into groups of a similar 
attitude. The content of work sheets is the same for all members of a group, but it can be different for 
different groups, according to your specific needs. If you find it useful, don 't hesitate to change a group 
and/or to demand for a curriculum related to your own questions on the problem. 

There is a collection of references quoted in the instruction. They are to be read only to such an 
extend which you find relevant and useful to solve the problems posed in the work sheets and to answer 
your own questions. 



3 Simplifying: analogy between Stokesian hydrodynamics and 
electrostatics 

Augustc Rodin: "The most remote antiquity is my habi- 
tat. I want to link the past to the present; to return to 
memory, judge it, and contrive to complete it. Symbols 
are the guidelines of humanity. They are no lies. "«*]/ 



3.1 Reasoning by analogy 



Read an introduction from [McDermott | , p. 90 



Electrostatics is simpler than Stokesian hydrodynamics. Therefore developing the analogy and show- 
ing 'how its corresponding parts are alike' helps to understand the basic concepts and processes of the 
complex technique we are going to study. To put the emphasis on foundations in this section we assume 
that there is no external ambient flow, which will be added to the system in Sec. [i[ 

3.2 Units 



To allow for an easy comparison with Kim, Karill^j and [Jackson], our reference textbooks, we will use 



SI units in hydrodynamics and CGS units in electrostatics, i.e. we assume that k=l. 

How the unit of charge (so-called statcoulomb) is defined in CGS system? 
How does this unit relate to centimeter, second and gram? 
Calculate how many coulombs it is. 



Reference: / Jackson ], Appendix 4 



3.3 Basic equation 

The basic equation of both electrostatics and hydrodynamics can be written in a general form as: 

L * = s (1) 
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where ^ is a physical held to be found, Lq is a differential operator, and s is a known source distribution. 
In electrostatics \& is the scalar potential held <£> and s is the charge density p. 

Specify the operator Lq in electrostatics. Use CGS units. 

Specify the meaning of an unknown field '5 , a differential operator Lq and a given source s in Stokesian 
hydrodynamics. What are the similarities and the differences in comparison to electrostatics? 

Compare your analogy with analogies developed by other groups. 
The eq. needs to be supplemented by boundary conditions. 

Specify what do you understand as the fluid boundries in the friction problem. 
References: [Jacksonl, Sec. 1.7, | Kim, Karille^ , Sec. 1.2.3. 



3.4 Uniqueness theorem 
Guiding question 

In electrostatics the solution to the Poisson equation is determined uniquely by specifying on the bound- 
ary: 

-the normal component of the electrostatic held E (so-called Neumann condition) or 
-the potential $ (so-called Dirichlet condition). 

Predict what is a hydrodynamic analogue of this theorem. 
3.4.1 Electrostatics 



Use the Green's identity ([Jackson], Sec. 1 



f OV 2 ?/' + V<t> ■ Vi/>)d 3 r = <f cj)V4> ■ ndA (2) 
Jv Js 

(where S is a surface surrounding V, n is the unit vector normal to S) to show the uniq ueness of solutions 
to the Poisson equation (|l|) for Dirichlet and Neumann boundary conditions (| Jackson |, Sec. f .9). Specify 
what is the meaning of uniqueness in both cases. 



3.4.2 Stokesian hydrodynamics 

Develop a similar proof for the Stokesian hydrodynamics. 
Use the Gauss theorem for a tensor K: 



d l K l ...id i r = d> Ki.,.imdA (3) 
v Js 

to formulate a useful generalization of the Green's identity (^|) for vector functions "0, <fi. 

State the uniqueness theorem. Specify what do you mean by uniqueness. 

Hint: Find out a vector analog of the following scalar theorem used in electrostatics: 

If Vj> = 0, t hen $ = cons t(r). 

References: [ Kim, Karilldj , Exercise 2.1 and Sec. 2.2.1, [ Pozrikidis] , Sec. 1.5. 



Compare your findings with results obtained by other groups. 
Answer the guiding question. 

3.5 Boundary conditions 

In Stokesian hydrodynamics from now on we will restrict to the stick boundary conditions, i.e. to the 
Quid velocity v(r) at the boundary equal to the rigid motion velocity of the boundary itself: 

v(r) = U + nxr (4) 
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However, according to iFcldcrhof 1976b], Fcldcrhof 198t], the formalism is valid for a more general class of 
the so-called slip boundary conditions, namely 



t-v(r) = t- (U + 12 x r) +lt- <r(r) • n 
n ■ v(r) = n - U 



at the boundary S 



where t is a unit vector tangential to the boundary surface, I vary from (for the stick boundary conditions ) 
to oo (the so-called perfect slip boundary conditions), <y is the fluid stress tensor: cy = fi(diVj +djVi) —pSij. 

A supplementary problem: slip boundary c onditions In Stokesian hydrodynamics the slip 

boundary conditions on the sphere surface are defined as fpcldcrhof 1976b\ : 

It ■ cr ■ n (5) 

(6) 



i • v 
n ■ v 



t ■ (U + fi x r) 
n-U 



where t is a unit vector tangential to the boundary surface, I vary from (the so-called no slip or stick 
boundary conditions) to oo (the so-called perfect slip boundary conditions). 

Do slip conditions determine a unique solution to the Stokes equations? Prove your statement. 



In Sec. 3 we assume that the fluid is motionless at infinity: v\oo = 0. In the following sections we will 
consider the general case of any conditions at infinity 

Construct simple electrostatic analogues of: 

A) (vi, fli) - a pair consisting of translational and rotational velocities of a body i, 

B) a body at rest in a viscous fluid, 

C) a rigidly moving body in a viscous fluid, 

Specify analogous equations for the boundary conditions in all cases. 



3.6 Friction problem 

The friction problem for N bodies is the following: 

If given: Cl a , U a , v (r), then what are Fp, Tp (a,(3— 1,...,N)? 

In this section we assume that the ambient fluid flow vanishes: Vo(r) = 0. In such a case motion of N 
bodies with velocities Up + dp x (r — r^) result in forces F Q and torques T Q exerted on the fluid by the 
body a, determined by the N-particle friction matrix £: 




Explain how does ([7|) follow from Stokes equations. Does £ depend on: 

- position in the fluid, r, 

- position of a body (3, Vp, 

- translational velocity of a body /3, Up and 

- rotational velocity of a body (3, ftp? 

If yes, explain how. If no, why not? What is the dimension of C? 

Are the a(3 components of the 2-particle friction matrix equal to C Q( 3j t ne corresponding components of 
the N-particle friction matrix from eq. (^)? Support your answer by a reasoning. 

Develop an electrostatic analogue of the friction problem (make it as simple as possible). What are 
the electrostatic analogues of the quantities appearing in eq. (]?])? Explain. 

In Sec. [|-fi| we consider N rigid bodies of an arbitrary shape. Later we will concentrate on N spheres 
only. 



3.7 Green function 

Green function G will help us to solve our friction problem, in a similar way it helps to solve its electrostatic 
analogue. 



G 



Guiding question 



Predict if the following statement is always true, true only under special supplementary conditions (if 
yes, specify them) or false; explain your reasoning: 

*(r) = J d 3 r'G(r,r')s(r') (8) 

3.7.1 Definition 

Green function G is a solution to the equation: 

L (r)G(r,r') = <5(r-r') (9) 

Is G a scalar, a vector or a tensor in: A) electrostatics; B) hydrodynamics? 
Write down (||) explicitely for electrostatics and for hydrodynamics, indicating arguments and all com- 
ponents. 



Specify G in electrostatics ([ [Jackson , Sec. 1.10) and in hydrodynamics ([Kim, Karillaj, Sec. 2.4.1) if 



the Dirichlet boundary conditions vanish at infinity. Make your definition consistent with your choice of 
From now on we will assume that G is the Green function for an infinite system. 



However, the formalism presented in Sec. H-W has been develop ed for any G reen function G \Fcldcrhof 198t 



corresponding to a container or to periodic boundary conditions { Hasimotc J, \Fcldcrhof 198$ , \Cichocki.Fcldcrhof 1989b) . 



3.7.2 Derivatives of the Green functions Q and "P 

The Green functions Q and "P satisfy the following identities: 

Liv 2 g t] (R) - diPj(R) = -M 3 (#) (io) 

diQij{R) = (11) 
V a Vj{R) = djS 3 (R) (12) 

where R — r' — r and all derivatives are taken with respect to r' : d L = djdr\ 



djPi(R) = S 3 (R). (13) 



Although Q and "P are functions, but their derivatives are distributions. We need to have a clear 
prescription how to evaluate such derivatives. Each of them can be understood as a limit of a sequence 
of functions. 



Show that the Green functions Q and V for the infinite system (our choice from Sec. 3.7) 



g l3 (R) = — f'% • ^-tM ii4i 



8irii \R R 3 

V * R) = 47§ + ^' (15) 
can be obtained as the following limits: Q = linia^o Q a and "P = lim a _ > o'P ) where 

QU R ) = ^(-9id j + 5 ij V 2 )(R 2 + a 2 ) 1 / 2 (16) 

VUR) = ~djV 2 (R 2 + a 2 ) 1 / 2 (17) 

Explain the procedure how to evaluate DQ and DT 3 for a differential operator D. 
Apply this prescription to verify eq. (11). 



3.8 Boundary integral equations 

Boundary integral equations are useful if there is a closed surface S surrounding a volume V and the held 
VP is defined both inside and outside. 
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3.8.1 Electrostatics 



Use the Green's theorem: 



(0V> - tpy 4>)d 6 r -- 
v JS 



dn dn 



dA 



to derive the following expressions for the electrostatic field potential $ ( | Jackson |, Sec. 1.10) 



/ p(r')G(r,r')d 3 r' + -j- I 

JV 47T Jg 



On' 



dn' 



dA' 



$(r) r eV 
r £ V 



Reference: [Jackson], Sec. 1.6, 1.8 



What is the hydrodynamic analogue of the electrostatic -^p-? 
Does (|l9| ) answer the guiding question? Explain. 

3.8.2 Stokesian hydrodynamics 

If there is no ambient Row, the hydrodynamic analogue of identities is given as: 



(18) 



(19) 



X 



pfjG 3 kd 3 r' 



{Gjkoij - Vjl^dlGjk + d'jGik) - SjiVk^n^dA' 







p(r) 




r e V 
reV 



(20) 



(21) 



Are <r and v taken at r or r'l Put the missing order of arguments of the Green functions Q and "P 
- (r, r 1 ) or (r 1 , r) - into eqs 

What are the symmetry properties of the Green functions Q and "P? Which of them are general, and 
which are due to the specific symmetries (no fluid motion at infinity) of the Oseen functions Jl4|)-(15)? 

How does the unit vector n! point: out or into the fluid? 

How do the equations simplify if there is no external forces acting on the fluid other than gravity? 



Reference: j Pozrikidis f, Sec. 2.3, Kim, Karilh ], Sec. 2.4.2, , Happel,Brcnnci J, Sec. 3.4. 

Eqs ©-© are valid if there is a closed boundary of any shape inside a fluid. How would you 
modify them to describe a rigid body in a fluid? Explain. 

The integral representation (j2^)-(j2l|) still does not allow to address the quiding question from Sec. 3.1 
-just the opposite, it seems to contain a term differing in form from eq. (0). We will come back to this 



problem in Sec. 4.2, but first we will learn how to take into account the existence of an ambient How 
around a particle. 



4 Particle in ambient flow 



Augustc Rodin: "First, I usually create my stone children 
without cloths. Then all I have to do is to throw some 
drapery over them...^ " 



4.1 Ambient flow 

How to use the results of Sec. || to construct solutions of the Stokes equation in the presence of an ambient 
flow? 
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4.1.1 Definitions 



1. The ambient How is a solution of the Stokes equation with given boundary conditions. 

2. The external ambient How Vo(r) is a solution of the Stokes equation with a given boundary condition 
at inHnity: 

vo(r)\oc=V{O,0) (22) 
Specify what is the ambient flow for the boundary condition: 

t>o(r)|oo =Vq. (23) 

Explain your reasoning. 

Give examples of other external ambient flows. 

Construct example of an ambient flow which has the same boundary conditions at infinity as a certain 
external ambient flow but which is different. Predict how this construction can be in agreement with the 



uniqueness theorem from Sec. 3.4. (We will come back to this example in Sec. 4.2.1 



4.1.2 Equvalence of solutions 

Assume that v(r) is the solution of the Stokes equation with the boundary conditions at infinity and on 
closed surfaces S a , a — 1, TV: 

v(r)\ Sa = U + nx(r\ Sa -r a ) (24) 
wtoloo - V(6,d>) (25) 

It is often said that the bodies S a are immersed in the ambient How Vo{r) given by (|2. 

Does v(r) — Vo(r) satisfy the Stokes equation? If yes, then specify the boundary conditions and 
explain your reasoning. If no, then why not. Explain how two solutions - in the presence and in the 
absence of an ambient flow - are equivalent . 

4.2 Formalism of induced forces 



Generalize the eqs (p0|)-(f2~ij) from Sec. 3.8.2 to describe a closed boundary in an ambient external flow 
vo(r). Explain your reasoning. 



The starting point of j Felderhof 1988\ (eq. 2.7) is the conjecture that the formal solution to the Stokes 



equation can be written as: 



v(r)-v (r) = J g(r,r')t(r')d 3 r' (26) 

where f(r) is the total force density exerted on the Huid. 

To understand and to justify this statement we will derive a formula in the form of (|2q) from the 



boundary integral equations given in Sec. 3.8.2 . To this goal we will Hrst develop the concept of induced 
surface force density. 

4.2.1 The concept of induced forces 

Construct a hydrodynamic analogue of the electrostatic relation between the induced surface charge 
density and the electrostatic field. How do you interpret the meaning of the adjective induced describing 
a surface density on the boundary in both cases? 



The goal is to express the boundary integrals given in Sec. 3.8.2 in terms of a surface induced force density 
(depending on the Huid stress tensor at the surface) rather than in terms of the Huid velocity v at the 
boundary. 

The question is how to achieve it - the eqs (^)-(^Jj depend on both a and v at the boundary. The 
idea is to Hrst introduce an artiHcial Huid How inside the rigid solid particles. Namely, to construct the 
inside solution to the Stokes equation with the same stick boundary conditions ( p^ j at the particle surface 
as those which determine the real outside solution. Next, to determine what are these induced forces, 



using the eqs ()2(\)-(W\), which you have just generalized for a non-vanishing ambient How. 



9 



Reference: [CoxJ3renner' ', Sec. 2 and , Mazur,Bcdeaux ], Sec. 3-4. 



Solve the Stokes equation inside a volume V with the stick boundary conditions ( |24| ) and the ambient 
flow given by (|25|). Hint: Make use of the uniqueness theorem. 

A supplementary problem: uniqueness theorem revised. The combination of the inside and the 
outside solutions is different than the ambient flow, although both satisfy the same boundary conditions at infinity. 



Explain how this construction can be in agreement with the uniqueness theorem from Sec. 3.4. Compare with your 
reasoning from Sec. |4.1.l| , where you also constructed example of an ambient flow which had the same boundary 
conditions at infinity as a certain different external ambient flow. 



4.2.2 Justification of the formalism 

Write down two sets of the generalized eqs (pp|)-(pT|): 

- for V being the interior of all the particles and 

- for V being the real fluid. 

Combine both sets and make use of the stick boundary conditions to eliminate the surface integrals 
including values of the fluid velocity at the boundary. 

Derive an expression for the induced force density in terms of the fluid stress tensor at the surface. 

How will you calculate the total force acting on a body in terms of the fluid stress tensor, and how 
in terms of the induced force density? 

How does the induced force density relate to f in eq. (|2^)? What is the range of integration in 
eq- ©? 

Use the generalized eqs (|2p|)-([2l|) to derive an equation for p(r) similar to eq. (prf). 

Specify the properties of the Green function used in this section. How could they be justified? 



4.3 Method of reflections 
4.3.1 The difficulty 

Assume that v„ is the solution of the one-particle friction problem, i.e. the solution to the Stokes equation 
with the boundary conditions ( p5| ) and the one-particle version of (pi|). 

Does v = y]g—i v Q satisfy Stokes equations? Explain your reasoning. 
Does v satisfy the boundary conditions ( |25| ) at infinity? If yes, then why? If not, then how could you 
construct v', another combination of v Q , satisfying (25)? 

Calculate v (and v') at the surface S a of the body a. How do they compare to the boundary conditions 
Are v (and v') solutions to the N-particle friction problem? Explain your reasoning. 



4.3.2 Construction 

Method of reflections is an iteration procedure to construct an approximate solution to the N-particle 
friction problem, building it from N single-particle solutions. At each step corrections are added to 
decrease descrepency between the boundary conditions and the actual value of the approximate solution 
on the body surfaces. It means that at each step we modify the single particle solutions by a better 
adjustement of their boundary conditions. The N-particle solution v is formally written as: 

N 

v = ^v a (27) 

a=l 

Each single particle solution v a for particle a is given by the following formal expansion: 

(28) 
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Each HJ 7 ,,, a satisfies the Stokes equations. 



Reference: j Kim, Karilla l, Sec. 8.1 



Specify what are the boundary conditions at infinity and at the surface S a for Vq, for w a , for wp a 
and for W.__fi a . 

What is the approximate value of v at the surface S a after: 

- the first 

- the second 

- the n-th 
iteration step? 

What is the effective ambient flow in which the particle a is immersed before: 

- the first, 

- the second and 

- the n-th 
iteration step? 

4.3.3 Interpretation 



Eq. ( |2iSP can be interpreted as a summation over all incident and outgoing "waves" in multiple scattering 
(in a sequence of reflections). 

Specify what are the incident and the outgoing "waves" scattered by a particle (3 at the first, the 
second and the n-th iteration step. 

Predict if the iteration procedure is convergent. Give arguments. 

Part III 

Exploring the structure: a lecture. 

Auguste Rodin: "For the first time I saw separate 
pieces, arms, heads or feet; then I attempted the figure 
as a whole. Suddenly, I grasped what unity was. . . "J*^/ 

5 Application of the basic concepts 



Reference: [ |Felderhof 1988fl , Sec. 2. 



5.1 Reformulation of the friction problem 

Forces and torques exerted by the fluid on the surface S a of the particle a are given in terms of the fluid 
velocity v and pressure p as: 

F Q = - <j> cr -n a dA (29) 
T Q = — <j) (r - r a ) x (cr • n Q ) dA (30) 
where v, p satisfy the Stokes equations: 

pS7 2 v-Vp = (31) 
V-v = (32) 

with given boundary conditions: 

v(r) -> V(e,<t>) r^oo (33) 
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v(r) = w a (r) = U a + n a x(r~r a ) r e S a (34) 

The surface normal n in (|2^)-(^o|) points into the particle. 

To solve (E9h-(p4j) we will apply the tools developed so far, in the following way: 



1. We evalute the ambient flow (vq, po) as the solution to (J3 1|) -(32) with 



2. We use the formalism of induced forces (Sec. 4.2) and the concept of equivalent solution (Sec. 4.1) to 



replace (gSj)-^ by: 

F a = J f d 3 r (35) 

T a = J(r-r a ) xid 3 r (36) 

w a (r)-v {r) = J g(r,r')i(r')d 3 r' for \r - r a \ < a (37) 

with Q given b y (|l4|) , and with f - non- vanishing on the particle surfaces only - related to the fluid stress 
tensor cr as in [Fclderhof 1976a : 



f(r) = £f Q (r) (38) 

a 

f Q (r) = a ■ n a S(\r - r a \ - a) (39) 
Eq. ([37]) has beed obtained by substitution of ( |34j ) into (p6|). 

5.2 Induced forces in terms of the boundary conditions for the fluid velocities 

The reformulated friction problem means solving ( |37j ) for f. First we note that due to linearity of the 
Stokes equations f depends linearly on the boundary conditions of the equivalent solution, i.e. the rigid 
motion of the spheres minus the ambient flow: w(r) — vo(r), where 

w(r) = J2 w «(r) (40) 

a 

and w a are given by (|34|). We write it as: 

f(r)= [ Z(r,r')[w(r r ) - v (r')]d 3 r' (41) 



and Z vanishes if r or r' is located outside a particle. 

Z has a matrix form, with Z a /3 relating f Q , the forces acting on the particle a, to the boundary 
conditions wp — vq, on the surfaces of all the particles j3. Using a simplified notation we write it as: 

f a = Z a p[wp-v Q ] (42) 

Therefore our goal is to find the N-body friction kernels Z a p(r,r'). We will do it in two steps: 

1. Simplification: one particle a in the fluid flow. We will have a look how to evalute the one-particle 
friction kernel Zq(oi). 

2. Multiplication: many particles. We will see how to express Z in terms of one-particle friction 
kernels Zo(a). 

5.3 Multipole expansion 

Assume that we know the N-particle friction kernel Z a p. The question is how to find the forces and the 
torques. Substituting ( (4l"| ) into ( p5| ) and ( |3^ ) we get: 

F a - J d 3 r J d 3 r'Z al3 (r,r')[w p (r')~v (r')} (43) 

T Q = J d 3 r J d 3 r'(r - r Q ) x Z a p(r, r'){w p {r') - v (r')] (44) 
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If there is no ambient flow, then: 

w a (r) - v (r) = U a + fl a x (r - r c 

and we have: 



C c y Q/3 V n 



The friction matrix elements can be written as: 

< b ia \Z af 3\bj(3 > 

with 



< 6oo| = l^oa > = 1 a (r) 

<b la \ = \b la > = -e ijk (r - r a ) k 9 a (r) 



and the scalar product defined as: 



< a b >= J a*(r) ■ b(r)d 3 r 



(45) 



(46) 



(47) 



(48) 
(49) 



(50) 



If there is no ambient flow, then: 



Wp(r)-v (r) = Y^ 



V p [wp{rp) - v (r p )] 



(51) 



We construct a complete set of functions bjg - combinations of (r — rp) p and a complete set of "velocity 
multipoles" C P p - combinations of V p [wp(rp) — Vo(rp)]. Instead of eq. ( p6| ) we now have: 

F a = < b 0a \Z a p\b p p > C p (52) 
T a = < bi a \Z al3 \b pf 3 > Cp (53) 



or equvalently: 



In particular: 



T„ 



^rt ^rr 



a/3 



r p 



Cp = Up 



(54) 

(55) 
(56) 



The questions remain how to choose the basis bp, how to expand the N-particle friction kernel Z and 
how to truncate. First we need to get acquinted with the N-particle and the one-particle friction kernels 
Z and Zq. 



5.4 Single particle solution 



The single particle solution in an ambient flow Wq is given as (26), 



v(r) - w (r) = J J g(r,r")Z (r",r')d 3 r" 



[w(r') - w (r')]d 3 r' 



(57) 



Note that r',r" are inside the particle, while r has no such restriction. To shorten notation we write 



(41) and (E7J) as 



f = Z [w-w ] 
v-w = QZ [w-w ] 



(58) 
(59) 
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5.5 Multiple scattering 



We can describe our many-particle system as the system consisting of a sphere a in an ambient flow 
created by the other bodies. But this ambient flow is also unknown, and it depends on the position and 
the motion of the sphere a. To solve this problem we will construct an iteration procedure. At each step 
the velocity field will be approximated by the sum of all one-particle solutions v a in a given ambient 
flow, evaluated from (|5^). Since this is not a many-particle solution, then vp will change the ambient 
flow in which particle a is immersed. It will be taken into account through modification of the ambient 
flow entering the next step of the iteration. 

The interpretation of this iteration in terms of a "multiple scattering" (or "reflection" ) on each particle 



f3 has been made in Sec. 4.3.3 



To carry out the multiple scattering we need to specify what is the total ambient flow w n , in which 
the particle a is immersed after each step n of the iteration procedure. w n is the "wave" outgoing from 
the step n and incident to step n + 1. It consists of the incident ambient flow w n -\ and corrections 
coming from the other particles (3 ^ a, evaluated from eq. (|59|) in step n — 1. Therefore: 



w 



W, 



W 7 



(60) 



- w n _i, n — 1, 2, ... correspond to the subsequent terms in the scattering expansion given in eq. (p8| 



In Sec. [4. 3.2 we have already analyzed the boundary conditions for them. Therefore, with the help of 



( 59 )-(p0[) the multiple scattering process made in step n can be described as: 

Vla-Vo = QZ 0a [w a -V ] 
V( n +l)a ~W n = -QZ 0a (w n - W n -l) 

The multiple scattering (|60| ) -(|6^) is equivalen t to the construction of the fluid velocity v by the metod 
of reflections [Kim, Karillaj, Happcl, Brenner |. Eqs (|60|) , (J6l]) and (p3) specify the ambient flow: 



for n > 2 



(61) 
(62) 



Wx-Vq 



If 



71+1 



- ^ G Z 0pi w n - Wn-l) 



for ri > 2 



and eq. 



8|) gives the induced forces: 

fla = Z 0a [w a — Vq] 



Z 0a (w a -w n ) = i na + Z 0a (lV n -i - w n ) 



for n > 2 



(63) 
(64) 



(65) 
(66) 



I (n+l)a 

Therefore the multiple scattering expansion of f a reads to the following form of the N-particle friction 
kernel Z a p in eq. (fl2"|): 



Z a p — Z 0a S a p — (1 — S a p)Z 0a QZ 



0/3 



qz q ~ ( qz 



Off 



(67) 



Auguste Rodin; "Yes, form I have looked at and under- 
stood, it can he. learnt: but the genius of form has yet to 
be studied. "1*1/ 
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Appendix. Announcement about the course 

Between June 15 and July 10, 1998 I will guide a course on 

HYDRODYNAMIC INTERACTIONS BETWEEN MANY SPHERES 

4 sessions, 2.5 hours each. 

GOAL: 

To inquire the basic structure and tools of the modern theory, which has been developed by Felderhof, Jones, 
Cichocki, Schmitz and their coworkers for 20 years, and resulted in numerical packages allowing for accurate 
calculations of hydrodynamic interactions. 

SUBJECT: 

How to determine the behavior of N spheres in low Reynolds number incompressible flow (N between several and 
several hundred). Namely: 

If an ambient fluid flow and external forces & torques acting on them are given, then what are their translational 
& rotational velocities (mobility problem)? 

If their translational & rotational velocities and an ambient fluid flow are given, then what are the forces & torques 
they exert on the fluid (friction problem)? 

KEY WORDS: 

Stokes equations, stick boundary conditions, ambient flow, Green function, induced forces, friction kernel, gen- 
eralized resistance matrix, generalized mobility matrix, multiple scattering expansion, force multipole moments, 
vector harmonics, rotational invariance. 

IDEA: 

To simplify, but knowing how to reach for complexity. 

Since the course is meant to be a first step needed to be done before making a more sophisticated analysis, then 
we will concentrate on basic concepts applied to a simple system. In particular, the following problems treated 
by this theory will be mentioned, but will not be discussed: slip boundary conditions, Green function other than 
Oseen tensor, lubrication phenomena, averaging procedure leading to evaluation of transport coefficients, mobility 
and friction problem for non-spherical shapes of particles (built from spheres). Analogy with electrostatics and 
quantum mechanics will be outlined, since you may later find it helpful in earring out calculations. 

ATTITUDE: 
To make it useful. 

Therefore first of all you are welcome to participate in making a plan of the course, by e-mailing me your sug- 
gestions what you would like to gain, what do you need it for and which concepts from those listed above are of 
your interest and which are not. Secondly, to help you in applying the technique, the sessions will be based on 
your active inquiry in small groups and your own solving of some basic problems rather than on passive listening 
to a lecture. Finally, your comments on time allocated to this activity are appreciated. 

CONTACT: 

If you want to participate, reply by e-mail before Thursday, June 11. Please let me know what are your time 
limitations - it will help me to fix the day of the week and the hour of our sessions. 

Maria Ekiel-Jezewska 
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